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ABSTRACT 

We quantify the detectability of stellar Milky Way satellites in the Sloan Digital Sky Survey (SDSS) 
Data Release 5. We show that the effective search volumes for the recently discovered SDSS-satellites 
depend strongly on their luminosity, with their maximum distance, D max , substantially smaller than 
the Milky Way halo's virial radius. Calculating the maximum accessible volume, V ma x, for all faint 
detected satellites, allows the calculation of the luminosity function for Milky Way satellite galaxies, 
accounting quantitatively for their detectability. We find that the number density of satellite galaxies 
continues to rise towards low luminosities, but may flatten at My ~ —5; within the uncertainties, the 
luminosity function can be described by a single power law dN/dMy = 10 x io 01 ( Ml '+ 5 ) ; spanning 
luminosities from My — — 2 all the way to the luminosity of the Large Magellanic Cloud. Comparing 
these results to several semi-analytic galaxy formation models, we find that their predictions differ 
significantly from the data: either the shape of the luminosity function, or the surface brightness 
distributions of the models, do not match. 

Subject headings: Galaxy: halo - Galaxy: structure - Galaxy: formation - Local Group 



1. INTRODUCTION 

In Cold Dark Matter (CDM) models, large spiral 
galaxies like the Milky Way and M31 form within ex- 
tensive dark matter halos from the merging and accre- 
tion of smaller systems. Although CDM models have 
had many successes on larger scales, one of the most seri- 
ous challenges facing CDM models is the so-called "miss- 
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Fig. 1. — The number of stars brighter than r ~ 22.5 in random 
realizations of Milky Way satellites of luminosity M r ~ —3, —5, —7 
(from bottom to top) with M92-like stellar populations, as a func- 
tion of distance from the Sun. The approximate number of stars 
required for a significant detection (by the algorithm described in 
Section |2) is ~30. 
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ing sa tellite " problem. First id entified by iKlvpin et al.1 
(| 1999T ) and IMoore et al.l (|1999f ). the problem manifests 
itself through the prediction by CDM models of at least 
1-2 orders of magnitude more low-mass sub-halos at the 
present epoch compared to the observed abundance of 
dwarf galaxies surrounding the Milky Way and M31. 

There have been a number of theoretical proposals to 
solve this problem. For example, the satellites that are 
observed could be em bedded only in the r arer, more mas- 
sive dark sub-halos (jStoehr et al.ll2002T ). or, the satel- 
lites may form only in the rare pe aks of halos that were 
above a given mas s at reionization (|Diemand et al.ll2005t 
IMoore et al.l l2006h ■ Alternatively, star formation in low 
mass systems coul d be inhibited by photoionization in 
the early Uni verse ([Bul lock et al.l l200ll ISomervilldl200l 
iBenson et al.l [2002) . All these ideas do not alter the 
abundance of dark matter sub-halos, but propose to solve 
the observed discrepancy by producing a smaller number 
of directly observable satellites, thus breaking any simple 
relationship between mass and luminosity. 

The known Milky Way dwarf spheroidal (dSph) satel- 
lites have been discovered by a variety of methods. The 
first seven were discovered serendipitously by visual in- 
spection of photographic plates, the Sextans dSph was 
found using automated scans of photographic plates and 
the Sagittarius dSph in a radial velocity survey of the 
Milky Way bulge. All-sky photographic surveys cover 
most of the sky away from the Zone of Avoidance, but 
searches of plates are limited to surface b rightnesses of 
~ 25.5 mag arcsec -2 (|Whiting et al.ll2007l ). The sample 
of known dSphs has long been bedeviled with selection 
effects, which are difficult to model with any accuracy. 
This situation has changed recently with the advent of 
very large area, homogeneous, photometric surveys such 
as the Sloan Digital Sky Survey (SDSS; York et al. 2000). 
The SDSS makes it possible to carry out a systematic sur- 
vey for satellite galaxies, which are detectable through 
their resolved stellar populations down to extremely low 
surface brightnesses. In essence, SDSS greatly facilitates 
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Fig. 2. — Differential convolution kernel applied to the stel- 
lar catalog to identify overdensities of a particular scale. A one- 
dimensional slice of the two-dimensional kernel is shown, where the 
width, or a, of the inner Gaussian is 6' and of the outer Gaussian 
is 60'. 

systematic searches for overdensities of stars in position- 
c qlor-magnitude s pace. 

IWillman et~aT1 (|2002j ) carried out the first SDSS- 
based survey for resolved Milky Way satellites, subse- 
quently discovering a new dwarf galaxy Ursa Major 
(|Willman et al.ll2005al) as well as an unusually lar ge glob- 
ular cluster, Willman 1 (jWillman et alj l2005bt) - al- 
though later evidence may favor its interpretation as a 
dark matter do minated dwarf galaxy with multiple stel- 
lar populations dMartin et alj l2007f). The color image 
"Field of Streams" (|Belokurov et all l2006a[ ). composed 
of magnitude slices of the stellar density in the SDSS 
around the North Galactic Cap, proved to be a treasure- 
trove for dwarf g alaxies, as Canes Ven at ici, Bootes I and 
Ursa M ajor II (|Zucker et all l2006allbl ; IBelokurov et all 
2006b) were all found in quick succession. A system- 
atic search in the "Field of Streams" led to the discovery 
of five more satellite galaxies, Canes Venatici II, Leo IV, 
Hercules, Coma, and Leo T, as well as another large glob- 
ular cluster, Segue 1 (IBelokurov et al.ll2007t llrwin et ail 
120071) . Very recently. IWalsh et alj (|2007l ) discovered an- 
other low luminosity satellite, Bootes II. 

As the faintest Milky Way satellites currently consti- 
tute our best markers of sub-halos, the faint end of the 
satellite luminosity function of the Milky Way satellites 
can provide stringent constraints on the process of galaxy 
formation, and can distinguish between a number of dark 
matter, structure formation and reionization models. So, 
it is important not merely to carry out a systematic sur- 
vey of the star overdensities in SDSS data for the dis- 
coveries per se, but also to compute the detection limits. 
These detection limits are the basis for a volume cor- 
rected luminosity function estimate and ultimately for a 
quantitative connection of satellite frequency and sub- 
halo abundance. Such is the purpose of this paper. It is 
important to note that for a volume-corrected estimate 
of the luminosity function, it is not nec essary to use ex- 
actly the same detection algorith ms as IBelokurov et al.1 
(j2007f ) or IWillman et all (|2005af ). Similarly, the detec- 
tion scheme does not need to be optimal for every indi- 
vidual dwarf galaxy. Provided the automated algorithm 
is able to detect all the Milky Way satellites, and the 



completeness properties of the algorithm are quantified, 
an estimate of the true luminosity function can be de- 
rived. 



2. DETECTION OF SATELLITE GALAXY CANDIDATES IN 
SDSS DR5 

The SDSS Data Release 5 (DR5) covers - 1/5 
of the sky, or ~ 8000 square degrees around the 
North Galactic Pole. SDSS imaging data are pro- 
du ced in five photome t ric bands, g, r, i, and 

z (iFukueita et al.l Il996t iGunn et all Il998t iHogg et all 
l200lUAdelman-McCarthv et alJ^OOalGunn et all 2006) . 
The data are automatically processed through pipelines 
to measure photometric and astrometric properties 
dLupton. Gunn fc Szalavl 119991; IStqughton et aIT l2002l; 
iPier et alj l2003t iTucker et alj 1200(1 ! All magnitudes 
quoted in this paper have been corrected for red- 
de ning due to Galactic extin c tion using the maps 
of Schlcgc l" Finkbeiner. fc Davisl (|1998h . Sometimes it is 
convenient to report our results in the V band, for which 
we use t he transformation V = g — 0.55(<? — r) — 0.03 
given bv lSmith et~aT1 (|2002t ). 

The SDSS data with the source catalogs used in this 
paper was downloaded from the SAI CAS V irtual Ob- 
servatory data center 4 (|Koposov et al. I [2007) and was 
stored locally in the PostgreSQL database. To per- 
form queries rapidly on the large dataset, we used the 
Q3C plugin for the spatial queries (jKoposov fc Bartunovl 
l200l . 

All the recent SDSS discoveries of dSph around the 
Milky Way, bar Leo T, are not directly visible in the flux- 
limited images, but were detected as overdensities of re- 
solved stars within certain magnitude and color ranges. 
This makes it straightforward to automate a detection 
method and assess its efficiency. The essence of any de- 
tection algorithm is to count the number of stars in a 
certain (angular) region on the sky, satisfying specified 
color and magnitude criteria, and compare the number 
to the background value. The excess of stars depends on 
the satellite's luminosity and distance. For a given lu- 
minosity, the distance fixes the number of stars brighter 
than the SDSS limiting magnitude, which is given by an 
integral over the stellar luminosity function. A simple il- 
lustration of the detectability of objects with a luminos- 
ity function like that of M92 is shown in Figure [TJ The 
curves show the number of stars brighter than r = 22.5 
for satellites of three different absolute magnitudes. The 
maximal distance probed by surveys like SDSS is con- 
trolled by the apparent magnitude of the brightest stars 
in the satellite. For intrinsically luminous objects, like 
CVn I (My = —7.9) , we can detect stars at the tip of 
the red giant branch at distances of up to ~1 Mpc. How- 
ever, for satellites with many fewer stars, like Hercules 
(My = —5.7), the giant branch tip is simply not popu- 
lated and we can only detect objects at distances up to 
-300 kpc. 

To identify the excess number of stars associated with 
a satellite, a common approach is to convolve the spa- 
tial distribution of the data with window functions or 
filters 5 . To estimate the star density on different scales, 



4 http://vo.astronet.ru 

5 This idea has a long history, particularly in algorithms for 
searching for features and clusters in imaging data. Widely used in 
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we use a Gaussian of width a, that is, 

L{x, y, a) = I(x, y) * g(x, y, a), 

where 

1 / x 2 + y 2 
g{x, y, cr = - — - exp — r - 

and I(x, y) is the distribution of sources 

i{x,y) = ^2 5 ( x - x *-.y - Vi) 



(1) 
(2) 

(3) 




This allows us to see the stellar density distribution 
at different spatial scales. For example, structures with 
a characteristic size of 1' will be more prominent when 
the stellar map is convolved with a 1' kernel, and less 
prominent when the map is convolved with 10' and 0.1' 
kernels. The resulting "blobs", or overdensities, can be 
easily identified on the differential image maps, namely 



AL = L(x, y,ai)- L(x, y, a 2 ) 

= I(x,y) * (g(x,y,a 1 ) - g(x,y,a 2 )) 



(4) 



Such differential image maps are generally convolutions 
of the original distribution with the kernel, which is a 
difference of two Gaussians. A one-dimensional slice of a 
kernel is shown in Figure [H When we convolve the map 
I(x, y) with such a kernel, we obtain an estimate of the 
local density minus an estimate of the local background 
(L{x, y, 02)). This interpretation allows us to quantify 
the significance as 



S{x,y,a 1 ,a 2 ) = - 



AL 

Oh 



(5) 



where a\ is the variance of L(x, y, <Ji). 

o\ = Variance(L(x, y, <J\)) — 

= Variance(I(x, y) * g(x, y, cti)) = I(x, y) * g 2 (x, y, a x 
= ^2 I ( x uVj) 9 2 (x - x u y - yj,ai) » 
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Under the assumption that a 2 > > o\ and a Poisson dis- 
tribution of the initial set of datapoints, the variance of 
S(x, y) is unity. This fact allows us to use the map of 
S(x, y) to identify overdensities above a specified signifi- 
cance threshold. 

3. APPLICATION TO SDSS DATA 

SDSS's morphological parameters (jLupton et al.ll200lh 
derived from the imaging data allow robust discrimina- 
tion between stars and galaxies down to r = 21.5. For 

astronomy are kernel-based density estimation methods, in which 
the density is obtained by convolving all the data points (inter- 
preted as delta-functions) with smoothly d ecaying kernels, which 
can be Gaussians (see e.g. Silverman 1986). A variant of this is 
used for feature detection in digital images in so-cal led scale-space 
science jLindcnbcrg 1993, 1998; Baba ud et aHIT986T) . 



Fig. 3.— The segmentation of the DR5 area into 17 32° X 32° 
fields, used for the stellar overdensity search described in the text. 



2L5 < r < 22.5, the discrimination is still reasonably 
reliable, but it becomes increasingly untrustworthy be- 
low r = 22.5. Moreover, the ca talog is 95% complete at 
r = 22.2 (|Stoughton et al.ir2002D and drops quickly below 
this magnitude. At the faint end, the "stellar" catalog of 
unresolved sources is polluted by faint galaxies which are 
intrinsically strongly clustered. We will see shortly that 
the main task in providing a clean sample of dwarf galaxy 
candidates is removal of the extragalactic contaminants, 
for which we will employ the SDSS galaxy catalog. 

To proceed with the convolution (Eq. [4j, the DR5 
field of view is split into 17 segments as shown in Fig- 
ure [31 The division is for computational convenience and 
to minimize distortion in the gnomonic projections. In 
practice, we select stars and galaxies with a magnitude 
cut-off of r < 22.5. Due to the properties of the ker- 
nel, we expect edge effects at the boundaries of the DR5 
footprint and we discard all overdensities within 1° of a 
boundary. We use a color-cut of g — r < 1.2 and kernel 
sizes with o\ — 4' and a 2 — 60'. The color cut is chosen 
to be as conservative as possible as regards inclusion of 
= the tip of the red giant branch stars for metal-poor popu- 
lations, whilst the kernel size is of the order of the angu- 
lar size of the known dwarfs (see next section for details). 
The color magnitude cut used in this work may not be 
optimal for the detection of each individual dwarf galaxy 
(e.g. the isochrone masks should definitely work better), 
but the primary goal here is not to define an optimal 
(6) algorithm, but rather to develop a consistent algorithm 
that can detect known objects, for which the detection 
efficiency can be determined. Figure 2] provides an ex- 
ample of the application of the detection pipeline to the 
stellar and galaxy catalogs of SDSS DR5. The method 
successfully removes the varying background to leave un- 
derdensities (black regions) and overdensities (white re- 
gions). The sample field of view chosen for Figure [4] 
contains the already kno wn Milky Way satell ites Will- 
man 1 and Ursa Major I (jWillman et al.ll2005ai rbl). They 
are both recovered in the stellar map with significances 
of Sstar = 13.07 and 5.95 respectively. However, as we 
see in Figure 01 unresolved sources in rich galaxy clusters 
such as Abell 773 and 1000, visible as prominent over- 
densities in the galaxy map, also show up in the stellar 
map as significant peaks. 

In order to remove false positives caused by galaxy 
clustering, we need to understand the significance S && \ of 
overdensities in the map derived from the galaxy cata- 
log. Equation ([7]) does not hold, because the underlying 
distribution is no longer Poissonian (Figure [5]). The left 
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Fig. 4. — A sample 22° X 22° area in the convolved maps of the SDSS DR5 stellar (left) and galaxy (right) catalogs. The positions 
of objects Ursa Major I and Willman 1 are marked by circles. The positions of galaxy clusters Abell 773 and Abell 1000 are marked by 
diamonds, and demonstrate that galaxy clusters may lead to significant peaks in the stellar map. The linear diagonal structures seen in 
both images are caused by SDSS stripes. The images were produced using a kernel specified by <ri = 4' and <T2 = 60'. We reject peaks in 
the convolved stellar map if they coincide with significant peaks in the galaxy distribution. 
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Fig. 5. — Left: The distribution of pixel values in the convolved 
star map (solid line) and galaxy map (dashed line) for one of our 
17 fields in DR5. The Gaussian model curves with width of 1.0 
and 2.3 centered on zero are shown in red. Right: The standard 
deviation in the galaxy map normalized by a Poissonian standard 
deviation as a function of kernel size. 

panel shows the distributions of 5 s tar and <S ga i for all pix- 
els in the same field of view as Figure|U For the stars, the 
convolved source count distribution is almost a Gaussian 
with unit standard deviation, whilst the distribution for 
the galaxies is broader. The right panel shows how the 
width of the S g& i distribution grows with increasing ker- 
nel width as the convolution samples coherent structures 
on larger scales. To assign significance to the overdensi- 
ties in galaxies, we rescale S sa \, dividing by its standard 
deviation. 

Next, we remove obvious false positives by rejecting all 
objects within the region marked by dashed lines in Fig- 
ure [6l namely the intersection of the regions Sstar < 20 
and Sgai > 2. This removes most, but not all, the false 
positives caused by galaxy clusters. Additionally, there 
remains contamination from galaxies with large angu- 
lar size. The SDSS photometric pipeline mis-classifies 
HII regions and stellar clusters in these galaxies as stars. 
We remove the contaminants by cross-correlating with 
the positions of galaxies in the Third Reference Cata- 
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Fig. 6. — Distribution of Milky Way satellite detections in the 
•Sstar versus Sgai plane. The circles mark the known Milky Way 
satellites, the triangles are RC3 galaxies, and the rhomboids arc 
galaxy clusters. Objects towards the top left of the figure are likely 
the result of contamination by galaxy clusters or spatially extended 
galaxies. The decision boundary is shown as a dashed line; objects 
to the right and below the dashed line are selected as candidate 
satellites. 



logue (RC3) of Ide Vaucouleurs et all (|1991f h Even so, 
there still remain objects at a moderately high level of 
significance whose nature is unclear. Most of these are 
probably caused by galaxy clusters or photometry arti- 
facts, as judged from examination of Hess diagrams and 
SDSS image cut-outs, but there may still be a very small 
number of genuine Milky Way satellites. 

We detect all the known Milky Way satellites, except 
Boo II, in a catalog with magnitude limit r < 22.5, an- 
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Fig. 7. — Left to right: Hess diagrams for Candidates X, Y and Z listed in Table [T] 



alyzed using a kernel with <j\ = 4'. The most marginal 
detections are Leo IV and Ursa Major I, which have sig- 
nificances Sstar = 6.10 and 5.95 respectively. Objects 
above the threshold are listed in Table QJ and include 
three likely false positives, which are "Candidates" X, 
Y and Z. The Hess diagrams of these three detections 
are shown in Figure [7] The Hess diagrams offer little 
evidence to support identification of the candidates as 
genuine satellites. Deeper data are needed to provide 
definitive classification of the candidates but for the pur- 
pose of determining the satellite luminosity function we 
exclude all three candidates, as false positives, from fur- 
ther consideration. 

It is prudent to search for candidate satellites on the 
map convolved with different inner kernels, since the ker- 
nel biases the algorithm towards objects of a preferred 
size. Therefore, we performed a search on the map con- 
volved with kernels of 2' and 8'. In the former case, set- 
ting the significance to 5 sta r > 6.5 results in the detection 
of all objects except UMa I and no false positives; in the 
latter case, setting S^tar > 6.0 includes all objects except 
CVn II, Leo IV, L eoT, UMa I and no false positives. 

Boo II, found bv lWalsh et al.1 (|2007l ) . is problematic for 
our algorithm. Boo II contains a very sparsely populated 
giant branch, and so the brightest stars are sub-giants 
and turn-off stars at colors of g — r < 0.5. Given our 
preferred cuts, Boo II is undetected. It can nonetheless 
be found with our algorithm, but only by optimizing the 
color and magnitude cuts, for example, to g — r < 0.5 
and 21 < r < 23. 

4. APPLICATION TO SIMULATED DATA 

To test our detection algorithm, we carry out an ex- 
tensive set of simulations in which we add mock dwarf 
galaxy satellites and globular clusters to the SDSS DR5 
catalog. In particular, we add to the catalog the g and r 
magnitudes of stars from the simulated objects, at speci- 
fied right ascensions and declinations. These augmented 
catalogs are then fed through our automated pipeline, 
and the number of stellar overdensities with significance 
above the threshold is calculated as a function of dis- 
tance, size and luminosity. We explore how changes in 
the g — r color cuts and kernel sizes {a\ from Eq. [7J affect 
the efficiency of the algorithm. 

The g and r photometry of all simulated objects is 
based on that of the globular cluster M92. The left panel 



of Figure [5] shows the color-magnitude diagram (CMD) 
of M92, together with a ma i n-seq uence and red giant 
branch ridgeline from iCleml (|2006f ). to which we have 
added a horizontal branch ridgeline. From the r-band 
data, we construct a main-sequence and red giant branch 
luminosity function and approximate it with a smooth 
fit, as shown in the middle panel of Figure [5] We also 
determine the luminosity function for the stars on the 
horizontal branch ridgeline. We populate the ridgelines 
using the luminosity function. The choice is appropri- 
ate, as M92 (12 Gyrs, [Fe/H] » -2) is typical of the old, 
met al-poor populations i n the Milky Way satellites (see 
e.g. Ivan den Be rgh 2000). Additionally, we add a scat- 
ter in r- and g- magnitudes, derived from a fit to the 
errors in the SDSS point-spread function photometry, as 
illustrated in the right panel of Figure [8j 

The spatial distribution of stars in the simulated ob- 
jects is chosen to follow a Plummer law, which is a reason- 
able fit to most of the Milky Way dwarf spheroidals (se e 
e.g. Ilrwin fc Hatzidimitrioul 119951: iKlevna et al.ll2002f) . 
For ellipticities less than 0.5 - which corresponds to 
the most flattened of the SPSS discoveries, Hercules 
and U rsa Major II (Zu cker et al.ll2006bt iBelokurov et al.l 
I2007T ) - the detection efficiency of objects barely changes 
with ellipticity. The Plummer radius, luminosity and 
distance are chosen to cover uniformly in logarithmic 
space the following ranges: Plummer radius lpc < r/i < 
lkpc, luminosity —11 < M v < and heliocentric dis- 
tance lOkpc < D < IMpc. We generate 8000 galax- 
ies with random right ascension and declination within 
the DR5 footprint. We then split the simulated sample 
into 20 distance bins to eliminate overlap between sim- 
ulated objects. The stars from the simulated galaxies 
are added to the DR5 stellar catalog. Figure [9] shows 
mock CMDs for simulated objects matching the recently 
discovered dwarf galaxies Canes Venatici I, Hercules 
and Ursa Majo r II (Zucker et al, 2006a: Bclok urov et al.l 
[20071 : IZucker et al.ll2006bf ). These are good approxima- 
tions to the observed CMDs of these objects. 

In our simulations, we test several inner kernel sizes. 
The reason is that for a given distance, the kernel size 
gives rise to an optimum physical size of the detectable 
objects. For example, at a distance of 50kpc, a kernel 
size of 4' corresponds to a physical size of ~60pc. As we 
want our algorithm to be sensitive to objects of different 
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Fig. 8. — Left: M92 color-magnitude data from Clem (2006) used as a template for our simulated Milky Way satellites, together with the 
ridge line for the main sequence and red giant branch. The ridge line for the horizontal branch is our fit to Clem's (2006) data. Center: The 
observed luminosity function of main-sequence and red giant branch stars in M92, together with our model fit of the luminosity function 
used in the simulations. Right: The photometric errors of the SDSS r-band photometry and our model fit used in the simulations 



M v =-7.66 D=220kpc M v =-5.58 D=140kpc 



M v =-3.33 D=30kpc 




Fig. 9. — Simulated color-magnitude diagrams for hypothetical Milky Way satellite galaxies with proper ties close to those of Canes 
Venatici I, Hercules an d Ursa Major II — the actual color-magnitude diagrams of these galaxies are given in IZucker et all H2006al lbl) and 
IBelokurov et al.l [120071) . 



sizes, we use three different inner kernel sizes, namely 
a i = 2', 4' and 8'. An object is considered to be detected 
if it is above a threshold on the map convolved with at 
least one of the kernels (the threshold for the 2' kernel 
is 6.50, for the 4' kernel - 5.95 and for the 8' kernel - 
6.00, see Section [3]). We refer to this procedure as the 
combined kernel. This is equivalent to the algorithm used 
in the previous Section, because the list of detections for 
2', 4' and 8' kernels includes all the known dwarfs. 

Figure [TO] shows two-dimensional efficiency maps as a 
function of luminosity and size in seven distance bins 
spanning the range 8 kpc to 1 Mpc. For Figure [TO] we 
have used the color cut of g — r < 1.2 and the combined 
kernel, together with an outer kernel of size 02 = 60'. 
Black corresponds to zero detection efficiency, and white 
to unit efficiency. The locations of the known Milky Way 
globular clusters and dwarf galaxies in this parameter 
space are recorded as red triangles and circles. While a 
number of known objects lie well within the efficiency 
boundary, some of the recent discoveries lie close the 
boundary. It is evident that there is no steady gradi- 



ent in efficiency, but rather a steep boundary between 
detectability and non-detectability. In fact, the primary 
contribution to the finite-extent of the gradient visible in 
the Figures is produced by the significant extent of the 
individual distance bins (the width of the distance bins is 
0.3 dex). The pixel size in magnitude is 0.8, and in logrh, 
it is 0.3. This means that there are typically 10 objects 
in each bin and so we expect moderate fluctuations due 
to shot noise. 

As the efficiency changes so quickly near the boundary, 
and as several objects lie close to this zone, we carried 
out more detailed simulations on objects similar to the 
known dwarfs. We created 1000 Monte Carlo realiza- 
tions of each of the known dwarfs, and fed them into 
the pipeline. Table [2] lists the derived detection efficien- 
cies for each object. The detection efficiencies are > 50 
%, with the sole exception of Boo II, confirming our as- 
sertion that the known satellites possess high detection 
probabilities. 

For the regime in which objects are larger than the 
kernel size, some of the stars belonging to the satellite 
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Fig. 10. — Detection efficiency maps for Milky Way satellites, shown as a function of luminosity and size for different distance bins. 
White indicates 100% detection efficiency, bl ack indicates 0%. Red circles mark the locations of the known dwarf galaxies, red triangles 
the known globular clusters (data taken from Harris (1996)). Notice that many of the very recent SDSS satellite galaxy discoveries occur 
near the boundary, where the detection efficiency is changing rapidly. 



Koposov et al. 



TABLE 1 

Objects Detected and Their Significances 



it ascension 


Declination 


Sstar 


S al 
ga 


Name 


205.539 


28.382 


170.13 


8.24 


NGC 5272 


168.355 


22.148 


170.08 


19.16 


Leo II 


198.220 


18.159 


165.00 


16.48 


NGC5024 


152.100 


12.289 


123.86 


14.94 


Leo I 


199.104 


17.696 


122.09 


5.59 


NGC 5053 


211.359 


28.527 


121.91 


5.61 


NGC 5466 


229.006 


-0.130 


115.10 


19.48 


Pal5 


260.038 


57.914 


100.02 


15.64 


Draco 


250.426 


36.467 


94.10 


2.80 


NGC 6205 


322.483 


12.147 


87.58 


13.30 


NGC7078 


182.516 


18.544 


79.40 


7.58 


NGC 4147 


260.008 


57.765 


75.07 


11.27 


Draco 


114.534 


38.873 


72.28 


6.80 


NGC 2419 


323.212 


-0.865 


65.52 


-1.75 


NGC 7089 


187.670 


12.395 


59.69 


3.25 


NGC 4486 


202.011 


33.549 


44.12 


1.68 


CVn I 


187.419 


8.003 


37.18 


0.93 


NGC 4472 


149.834 


30.742 


28.13 


12.65 


Leo A 


190.698 


2.682 


27.73 


5.51 


NGC 4636 


114.608 


21.581 


25.88 


-6.25 


NGC 2420 


259.027 


43.063 


22.76 


-5.60 


NGC 6341 


183.904 


36.310 


20.43 


7.49 


NGC 4214 


185.036 


29.286 


17.28 


0.63 


NGC 4278 


210.010 


14.503 


16.95 


0.32 


Boo I 


190.773 


11.598 


16.84 


-1.52 


NGC 4647/4637/4t 


186.368 


12.909 


14.28 


0.76 


NGC 4374 


178.814 


23.371 


13.90 


19.34 


Abell 1413 


186.444 


33.539 


13.10 


16.06 


NGC 4395 


148.904 


69.081 


13.08 


19.85 


NGC 3031 


162.325 


51.051 


13.07 


0.11 


Willman 1 


242.741 


14.956 


12.52 


-0.07 


Pal 14 


186.315 


18.181 


11.59 


-3.72 


NGC 4382 


132.830 


63.124 


11.40 


0.92 


UMa II 


186.745 


23.913 


11.22 


-1.02 


Coma Berenices 


143.721 


17.058 


10.96 


3.85 


Leo T 


188.911 


12.544 


10.91 


-0.97 


NGC 4552 


210.691 


54.332 


10.67 


13.53 


NGC 5457 


151.369 


0.070 


10.64 


4.11 


Pal 3 


186.109 


7.294 


9.85 


-0.76 


NGC 4365 


172.319 


28.961 


9.53 


0.26 


Pal 4 


247.764 


12.789 


8.91 


1.05 


Hercules 


197.870 


-1.335 


8.23 


22.00 


Abell 1689 


194.292 


34.298 


7.39 


-4.61 


CVn II 


196.743 


46.569 


7.30 


8.29 


Abell 1682 


193.379 


46.415 


6.71 


-2.41 


Candidate X 


168.146 


43.440 


6.52 


-3.22 


Candidate Y 


352.182 


14.714 


6.39 


1.47 


Pegasus 


202.388 


58.404 


6.37 


-0.59 


NGC 5204 


225.323 


1.672 


6.34 


1.36 


NGC 5813 


187.038 


44.090 


6.13 


-2.49 


NGC 4449 


173.235 


-0.554 


6.10 


2.92 


Leo IV 


0.807 


16.097 


6.09 


-2.66 


NGC 7814 


179.144 


21.049 


6.05 


1.14 


Candidate Z 


184.843 


5.786 


6.04 


-1.25 


NGC 4261 


149.993 


5.316 


6.03 


2.12 


Sextans B 


139.470 


51.718 


5.97 


17.86 


Abell 773 


179.223 


23.379 


5.96 


5.22 


galaxy cluster 


158.695 


51.918 


5.95 


3.72 


UMa I 



-10 



-6 



Mv 



Fig. 11. — Characterizing the satellite galaxy detectability: 
this illustrative figure shows the model function e{My,^) (from 
Eq. [5} used to fit the observed detection efficiencies from the 
simulations and demonstrating the role played by the thresholds 
M\s lim an d Mlim- The function parameters used to produce the plot 
;8were Myjim = ~2mag, in i m — 29.5 mag arcsec ^, um — 1 mag, 
ct m = 1 mag 

surface brightness and absolute magnitude by adopting 
a functional form: 



e(M v ,n) = G 



/M 



v 



Mi 



V.lim 



CM 



G 



(8) 



where G denotes the Gaussian integral, which is defined 

as 



G(x) 



1 



exp(-r/2)dt. = -er/c 



(9) 



are missed by the window function, and for such objects 
the detectability is determined by the number of stars 
within the window function, i.e. the surface brightness. 
This effect produces the surface brightness limit seen in 
Figure [TO] For objects smaller than the kernel size, all 
the stars are within the window function regardless of 
the size of the objects, therefore for such objects, the de- 
tectability doesn't depend on physical size, but depends 
only on the total number of stars, i.e. the luminosity. 
This effect produces the rapid change in detection effi- 
ciency at fixed absolute magnitude evident in Figure [TTJ1 
These two regimes can be modeled with thresholds in 



To describe the detectability in each distance bin, there 
are four parameters that are fitted - namely the detection 
thresholds in surface brightness /ii; m and absolute mag- 
nitude Mv,iimi together with their widths c M and <jm- 
As an illustrative example, the grey-scale map of the ef- 
ficiency function e(My,/i) from Eq. [5] is shown in Fig- 
ure [Til with dashed lines indicating the thresholds. Note 
the shape of the detection boundary, with the prominent 
"knee" , which corresponds to the boundary between the 
two detection regimes for objects of different sizes at fixed 
distance, as described above. 

The two key parameters for the detection pipeline are 
the inner kernel size o\ and the color cut applied to the 
source catalogs. The top two panels of Figure [12] show 
-the dependence of Myfim an d fJ>Hm on distance, when con- 
volved with the three different inner kernels. For a given 
kernel, the limiting magnitude declines roughly linearly 
with the logarithm of distance. Objects at the limiting 
magnitude have an apparent size that is smaller than the 
kernel size and their detection significance is reduced by 
the background contribution. Shrinking the kernel size 
removes some of the background and increases the signif- 
icance of fainter satellites. The dependence of the Mv,Um 
and jttiim on distance for the combined kernel is not plot- 
ted, because in the top left panel of Figure [T2] the com- 
bined kernel basically follow the dependence of 2' kernel 
and in the top right panel of Figure [T^] the combined ker- 
nel follow the dependence of 8' kernel. This is illustrated 
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Fig. 12. — The known satellites and globular clusters shown in two-dimensional plots of Galactocentric distance versus absolute magnitude 
(left) and surface brightness (right). Circles mark the locations of the known dwarf galaxies, triangles the globular clusters. The error 
bars show either o-jy/ or cr^ derived from our model fits (see Eq. [8j. The detectability of the objects depends on their location relative to 
the limiting absolute magnitude (left) and surface brightness (right) as a function of Galactocentric distance for each kernel size/color cut 
employed in the search. Upper panels: The three lines show the detection limits for different sizes of the inner Gaussian in the kernel (blue 
— 2', green - 4', red - 8'). Lower panels: The four lines show the detection limits for the different g — r color cuts employed (black — 0.2, 
blue - 0.4, green - 0.6, orange - 0.9, red —1.2) and fixed kernel size of 4'. 



in the top left panel of Figure [T2l where it is clear that a 
smaller kernel allows us to detect fainter objects. How- 
ever, as the top right panel shows, this is at the expense 
of satellite size. Larger kernels pick up more stars from 
extended objects and hence reach fainter surface bright- 
ness. When combining different kernels in the pipeline, 
the overall limits in surface brightness and absolute mag- 
nitude (2', 4', 8', see Section© can be approximated by 
the blue line in the top left panel of Figure [T2] and the 
red line in the top right panel of Figure [121 Smaller 
kernels allow the detection of galaxies that are low in ab- 
solute magnitude, and larger kernels allow the detection 
of galaxies that are fainter in surface brightness. It is 
also reassuring to see that the error bars <jm are of the 
same order as the difference in the limiting magnitude 



moving to a neighboring bin. 

We explore the effects of changing color cuts and report 
the results in the bottom two panels of Figure Q21 The 
color cut of g — r < 0.4 can improve slightly the magni- 
tude limit for nearby objects by selecting turn-off stars. 
This improvement deteriorates rapidly as we exhaust the 
supply of turn-off stars. At larger distances, red color 
cuts like g — r < 1.2 are more efficient at picking up gi- 
ant stars. The same effect explains the drop in /j,n m and 
Mv,\im at large distances for bluish color cuts. Our choice 
of color cut g — r < 1.2 is conservative, mostly eliminat- 
ing thin disk stars, and is, overall, the best-behaved and 
most robust. It also allows us to minimize the influence 
of metallicity and age changes in the stellar population 
of the satellites. 
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TABLE 2 
Detection Efficiencies of 
Simulated Objects 
Resembling Known 
Satellites. 



Object Efficiency 




Bootes 


1.00 


Draco 


1.00 


Leo I 


1.00 


Leo II 


1.00 


SEGUE 1 


1.00 


Canes Venatici I 


0.99 


Willman I 


0.99 


Coma 


0.97 


Koposov 1 


0.90 


Leo IV 


0.79 


Ursa Major II 


0.78 


Leo T 


0.76 


Hercules 


0.72 


Ursa Major I 


0.56 


Koposov 2 


0.48 


Canes Venatici II 


0.47 


Boo II 


0.20 


TABLE 3 




Limiting Satellite Absolute 


Magnitude and Surface 


Brightness as a Function of 


Distance 




Distance My u m 




kpc 


mag/D" 


11 0.6 


27.5 


22 0.4 


28.7 


45 -1.9 


29.6 


90 -3.4 


30.0 


180 -4.4 


29.9 


260 -5.9 


29.9 


720 -7.5 


29.6 




-12 -10 -8 -6 -4 -2 
M v , [mag] 



Fig. 13. — The accessible volume within the DR5 footprint for 
galaxies with different lumi nosi ties and surface brightnesses /xii m , 
fjj < 30 mag/ □" (see Figure [12] The volume limited by the virial 
radius (280 kpc) and within DR5 is shown by the dashed line. 



Fig. 14. — The luminosity functions of Milky Way satellite galax- 
ies within ~280 kpc (virial radius) inferred from our analysis un- 
der the assumption of two different radial distributions of satel- 
lites, NFW-like (solid black line) and isothermal (dashed black 
line). The calculation uses the satellite list and the volume cor- 
rection factor obtained with the pipeline using the cuts r < 22.5 
and Sstar > 5.95. The arrows on error bars indicate that there 
is only one galaxy in the particular bin, and so the Poisson er- 
ror is formally 100% . The theoretical prediction of Figure 1 of 
Benson ct al. ( 200^) is shown in a red line, and the prediction of 
Somervillel (2002'l for z re i on = 10 is shown as a blue line. Addition- 
ally, the luminosity function for the bright (My < —11) satellites 
of the Milky Way sampled over the wh ole sky together w ith the 
bright M31 satellites within 280 kpc from lMetz et al'l (120071 ) is plot- 
ted with filled small symbols (the list of plotted objects consists 
of Sgr, LMC, SMC, Scu, For, LeoII, Leol, M32, NGC 205, And 
I, NGC 147, And II, NGC 185, And VII, IC 10). The function 
dN/dM v = 10 X 10 01 ( M v+5) is s hown in grey. 

5. THE LUMINOSITY FUNCTION 

5.1. Analysis of the detection efficiency maps 

With an understanding of which satellite galaxies can 
be detected in SDSS DR5, together with our sample 
of actual detections, we can now estimate the luminos- 
ity function of faint Milky Way satellites. We start by 
re-examining the efficiency maps in Figure 1101 where 
the locations of the known Milky Way globular clusters 
and dwarf galaxies are overplotted as red triangles and 
circles respectively. We can conclude that within the 
DR5 footprint there are certainly no bright satellites (ei- 
ther globulars or galaxies) nearby (D < 32 kpc) that 
have elu ded discovery . Howe ver, the disrupting galaxy 
UMa II (|Zucker et al.1 l20~06b| ) provides a clue as to the 
likely locations of remnants. It is still possible that dis- 
rupted galaxies remain undiscovered nearby. They can 
lurk in the black portions of the uppermost two panels 
of Figure [H 

All that has survived in the inner Galaxy (8 < D < 16 
kpc) is a population of globular clusters, which oc- 
cupy a small region in the luminosity and size param- 
eter space. They are predominantly old globular clus- 
ters belonging to the bulge. Only the densest sur- 
vive against the disruptive effects of Galactic tides and 
shocking, which is illustrated by the apparent size bias. 
Notice that the datapoints lie well away from the de- 
tection boundary, suggesting that the sample is com- 
plete at least within 8 < D < 16. Moving outwards 
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(16 < D < 32 kpc), the globular cl usters belong to the 
halo and may have been accreted (jMackey &; Gilmorel 
l200l . Their size distribution is broader. Some of these 
objects are i n the process of disruption, such as Pal 5 and 
NGC 5466 (lOdenkirchen et al.ll200ll Dehncn et al.ll2004 
iBelokurov et al.ll2006cT Fcllhaue r et al.ll2007t ). The very 
faint and distan t glob ular clusters discovered recently by 
IKoposov et al.1 (|2007t ) are visible in the third panel of 
Figure [TO] (32 < D < 64 kpc) right on the border of 
detectability. Further such sparse globulars may remain 
undetected. 

Beyond 30 kpc, the dwarf spheroidals begin to appear. 
The long-known dwarfs such as Draco and Sculptor lie 
far from the boundary, in regions of the luminosity and 
size parameter space where the DR5 search efficiency is 
unity. However, all the recent SDSS discoveries, such as 
Canes Venatici I, Bootes and Hercules, lie close to the 
detection boundary , where the efficiency de clines rapidly 
from unity to zero. Bel okurov et alj (|2007l ) claimed that 
there is a paucity of objects with half-light radii between 
~ 40 pc and ~ 100 pc. Our calculations support the idea 
that the gap is real and not produced by selection effects. 
If there were objects with radii between ~ 40 pc and 
~ 100 pc, there is a broad range of parameter space in 
which they would have been found. 

Although most of the new detections lie in the gray 
areas of the plot, the empty white regions with unit ef- 
ficiency are telling us something important. There are 
swathes of the parameter space in which we would have 
detected objects if they existed. For example, there are 
very few bright objects (M v < —6). The absence of de- 
tections of bright objects does by itself provide a strong 
constraint on the luminosity function of Milky Way satel- 
lites. There also do not appear to be any analogues 
of the extended , lumi nous star clusters found in M31 
bv iHuxor etal] (|2005h . Although SDSS data may still 
contain evidence for further, hitherto unknown, dwarf 
galaxies, it is unlikely that their nature can be unambigu- 
ously established without substantial quantities of follow- 
up imaging. We emphasize that, since we never probe 
fainter than a certain surface brightness limit, an even 
larger population of very low surface brightness galaxies 
- which can not be detected with SDSS - may exist. 

5.2. Estimation of the Luminosity function 

Figure [13] shows the accessible volume for galaxies of 
different luminosities probed by our algorithm (which in 
practice is a function mostly of the luminosity) within the 
SDSS DR5 footprint. As the logarithm of distance scales 
roughly linearly with limiting magnitude (see Figure [T2|) . 
so does the logarithm of the accessible volume. Using 
this, and the fact that the SDSS survey covers ~ 1/5 of 
the sky, we can convert the set of known objects into a 
volume corrected luminosity function 6 . 

The observed luminosity function is constructed using 
all the well-established dwarf galaxies in DR5, namely 
Leo II, Draco, Leo I, CVn I, Boo I, Hercules, UMa II, 

6 The existing data on the globular cluster population indicate 
that at least some globular clusters have complicated metallicity, 
age distributions a nd kinematics and may in fact be stripped nuclei 
of dwarf galaxies l|Zinnecker et al.|[T988t IPiotto et ^120071) . The 
selection of such objects which are considered as dwarf galaxies 
is an additional source of uncertainty in any luminosity function 
determination. 



Com, CVn II, Leo T, UMa I, Leo IV as well as the pos- 
sible dwarf Willm an 1. Segue 1 is not used because it 
is not in DR5 (|Belokurov et al.|[2007h . and Boo II is 
not used because it is not detected with our adopted 
identification-pipeline parameters. All the satellites in- 
cluded in our calculation have a surface brightness of at 
least 30 mag arcsec -2 . To relate the observed number 
of satellite galaxies in our sample to the total number 
of satellites in the Milky Way halo, it is necessary to 
adopt an underlying radial distribution of satellite galax- 
ies (see Appendix). In a given magnitude interval, we 
know the observed number of satellites within V max (M v ) 
from Figure [13] together with their detection efficiencies 
from Table [5J If we assume a number density law rt(r) 
for the satellites, then its normalization at each mag- 
nitude interval can be fixed by integrating the density 
law out to Umax- The total number of satellites within 
280 kpc (the virial radius of the halo) is now the inte- 
gral of the density law out to this limit. Figure [14] shows 
the results of the calculation for two such density laws. 
The dashed line shows the luminosity function assum- 
ing the satellites are distributed in an isothermal sphere 
(namely, n(r) oc 1/r 2 ). The solid line shows the lumi- 
nosity function if the density fall-off is steeper at large 
radii (ri(r) cx 1/r 3 , like Navarro-Frenk- White profile, al- 
though to prevent the 1/r 3 profile from diverging in the 
MW center we use n(r) cx r~ 2 (r + r c ) _1 with the core 
radius r c = 10 kpc). Of course, the nature of some of 
the objects we have included in the dwarf galaxy lumi- 
nosity function is still uncertain - in par ticular, Will- 
man 1 may be a globular cluster, although iMartin et al.l 
(|2007f ) provide evidence for a metallicity spread and dark 
matter content. It is unclear whether Leo T should be 
included or excluded, as it is most likely a transition 
object with rather different properties from the other 
dwarf spheroidal galaxies in our sample. The error bars 
in Figure [14] are given by the square root of the num- 
ber of datapoints in the absolute magnitude interval di- 
vided by the volume correction factor. At the bright end, 
the error bars are large, since we have only two objects 
with M v < —9, namely Leo I (M v = —11.5) and Leo II 
(M v = —9.6). At the faint end, the error bars are also 
large because of the substantial volume correction factor. 
In Figure [14] we show the luminosity function for satel- 
lites within 280 kpc (a proxy for a Milky Way virial radius 
(|Klypin et al.l l2002t iBenson et all 120021 )). To define the 
bright end of the luminosity function, which cannot be 
reliably determined from our data since DR5 does not 
contain dwarfs brighter than My ~ —11, we have also 
included in Figurc fMl thc estimate of the luminosity func- 
tion (filled points) for the bright satellites of the Milky 
Way sampled over the full sky, toget her with th e brigh t 
M31 satellites within 280 kpc from iMetz et alj (|2007l ). 
In Figure HU we also overplot the power-law function 
dN/dM v = 10 x io 01 ( M v+5) ; w hich approximates the 
datapoints in the range of— 19 < My < —2 (with prob- 
ably some flattening at My ~ —4). The integration of 
this power-law gives approximately 45 dwarfs brighter 
than -5.0, and 85 dwarfs brighter than -2.0. 

There are a number of theoretical predictions of the 
luminosity function of the Lo c al G roup in the litera- 
ture. For example, ISomerv illc (2002) shows the results 
of semi- analytic galaxy formation calculations, includ- 
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ing the effects of supernova feedb ack and photoion iza- 
tion. The lu minosity function from ISomervilld (|2002T ) for 
Zreion — 10 (|Page et al.ll2007l ) are plotted with blue line 
in Figure [Mj Although the numbers of luminous satel- 
lites are in reasonable agreement with the data, the shape 
of the luminosity function is not. All Somerville's (2002) 
luminosity functions turn over at M v « —9 or brighter, 
depending on the epoch of rcionization, whereas the lu- 
minosity function derived in Figure [H] turns over fainter 
than My w -5, if at all. Therefore, Somerville's (2002) 
theoretical calculations overproduce Draco-like objects 
{My —10) by a factor of a few, and underproduce 
much fainter galaxies like Boo (My ~ —6) by almost an 
o rder of magnitude. 

iBenson et al.l (|2002f ) also provides calculations of the 
luminosity function of the Milky Way satellites, includ- 
ing the effects of tidal disruption as well as photoioniza- 
tion. They report the luminosity functions for dwarfs 
with a range of of different central surface brightness 
cuts, namely 18, 20, 22, 24 and 26 mag arcsec -2 , the last 
of which is plotted in Figure [H] in a red line. At first 
glance, the fit seems plausible, especially given the size 
of the error bars on the datapoints. The turn over in Ben- 
son et al.'s luminosity function is at My w —3 and the 
numbers of predicted satellites at faint magnitudes are 
also consistent given the uncertainties. However, Benson 
et al.'s model significantly underproduces the number of 
bright satellites. Additionally, Benson et al.'s satellites 
have a much higher central surface brightness - our SDSS 
survey corresponds to a surface brightness c ut of ~ 30 
mag arcsec -2 . Figure 2 of lBenson et a l. (2002) does show 
the luminosity function for all objects, irrespective of sur- 
face brightness. Although there has been a large change 
in the luminosity function on moving from a detection 
threshold of 22 to 26 mag arcsec -2 , there is only a small 
change on moving from 26 to oo mag arcsec -2 . 

6. CONCLUSIONS 

There have been persistent discrepancies between the 
observed numbers of Milky Way satellites and the pre- 
dictions from numerical simulations of galaxy formation 
for a number of years. Although here has been a caval- 
cade of discoveries of new Milky Way satellites using the 
SDSS over the last two years, a systematic search - with 
quantifiable detection limits and efficiencies - not been 
undertaken. In this paper we have presented a quantita- 
tive search methodology for Milky Way satellite galaxies 
in SDSS data and have used this method to compute de- 
tection efficiency maps, which ultimately allow the con- 
struction of the satellite galaxy luminosity function. 

In our method, the star count map is convolved with 
a family of kernels which are the difference of two Gaus- 
sians. Intuitively, this algorithm can be understood as 
constructing an estimate of the local stellar density mi- 
nus the background. By attaching a statistical signifi- 
cance to the overdensities in the convolved image, this 
enables us to construct a ranked list of candidates. Al- 
though this idea is simple enough, its practical applica- 
tion is hampered by the fact that the separation between 
stars and galaxies by the SDSS pipeline becomes unreli- 
able at magnitudes fainter than r ~ 22.5. The resulting 
false positives must be removed by cross-correlating with 
galaxy catalogs. The significance threshold of peaks in 
our survey is set by requiring the detection pipeline to 



produce a "clean" list of Milky Way satellites. 

To compute the detection efficiency, we create mock 
SDSS catalogs with stars from simulated dwarf galax- 
ies and use Monte Carlo methods to estimate recovery 
as a function of satellite galaxy parameters and helio- 
centric distance. There is a sharp boundary between 
detectability and non-detectability. The efficiency maps 
make clear that there are large domains in parameter 
space in which objects would have been detected had 
they existed. In particular, even at heliocentric distances 
as great as 1 Mpc, objects brighter than M v ~ —8 would 
have been detectable in SDSS. Similarly, populations of 
extended, luminous star clusters would have been found 
in SDSS, if they existed in the Milky Way. 

With the efficiency in hand, we can - for the first 
time - correct the observed luminosity function of the 
Milky Way satellites for selection effects and compute the 
true luminosity function. The number density of satel- 
lite galaxies continues to rise well below My ~ — 8 m ; 
depending on the radial distribution model assumed it 
may or may not flatten or turn over at My > —5. 
Overall, the luminosity function of all Milky Way satel- 
lites may be reasonably well described by a power-law, 
dN/dMy = 10 x iq0-1(m v +5) from My = _ 2 to -18. 
This power-law predicts <~ 45 satellites brighter than 
My — —5, and ^85 satellites brighter than My = —2. 
The normalization of the luminosity function is in rea- 
sonable agreement with the predictions of semi-analytic 
modeling of galaxy formation, but the shape is not. 
There also remains a discrepancy in the distribution of 
surface brightnesses of such objects, in the sense that the 
semi-analytic models underproduce dwarfs with a central 
surface brightness fainter than 26 mag arcsec -2 . 
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APPENDIX 

THE CALCULATION OF THE CORRECTION TO THE LUMINOSITY FUNCTION 

To calculate the luminosity function of Milky Way satellites within tlf = 280 kpc, we select all the satellites within 
DR5 which are interior to tlf, and construct the histogram of My of these objects. From the simulations, we know 
that not all objects are detected with 100% efficiency and the histogram h(Mv) is weighted with the object detection 
efficiencies. 

h(M v )=J2-S(M v ,M v ,i) 



where is the detection efficiency of i-th object, My,i its luminosity and S(My, Mv,i) = 1, if My and My,i are within 
one bin of the histogram, and otherwise. 

Figure [T2l shows how the maximal accessible distance depends on the galaxy luminosity (the r max (My) function). 
From this function, we can construct the maximal accessible volume within the DR5 footprint (which covers 1/5 of the 
sky) as a function of galaxy luminosity, namely V mayi (My) = 47r/3 /d_rs rj^ ax (My)(see Figure[l3]), where Jdr5 is the 
fraction of the sky covered by DR5 . Then we construct the incompleteness correction c(My), using the probability 
distribution of the satellites n(r). When the maximal accessible distance for a galaxy is greater than tlf, the correction 
is 1, if not it is equal to the ratio of number of satellites within r max (My) to the number of satellites within r^-p: 



r max (Mv) 

J n(r)r 2 dr 



c(My) = < 



(I 



if r max (M v ) < r LF 



t-lf 

J n(r)r 2 dr 
o 

1 if 1 max 

(My) > 7"lf 



Finally, the luminosity function is obtained by dividing the histogram of luminosities h(My) by the incompleteness 
correction c(My) 



REFERENCES 



Adelman-McCarthy, J. K., et al. 2006, ApJS, 162, 38 

Babaud J., A. P. Witkin, M. Baudin, and R. O. Duda. 1986, IEEE 

Trans. Pattern Anal. Mach. Intell., 8, 1, 26 
Belokurov, V., et al. 2006a, ApJ, 642, L137 
Belokurov, V., et al. 2006b, ApJ, 647, Llll 

Belokurov, V., Evans, N. W., Irwin, M. J., Hewett, P. O, & 

Wilkinson, M. I. 2006c, ApJ, 637, L29 
Belokurov, V., et al. 2007, ApJ, 654, 897 

Benson, A. J., Frenk, C. S., Lacey, C. C, Baugh, C. M., Cole, S. 

2002, MNRAS, 333, 177 
Bullock J.S., Kravtsov A.V., Weinberg D.H. 2001, ApJ, 548, 33 
Clem, J. L. 2006, Ph.D. Thesis, 

Dehnen W., Odenkirchen M., Grebel E. K., Rix H.-W. 2004, AJ, 
127, 2753 

Diemand J., Madau P., Moore B. 2005, MNRAS, 364, 367 
Fellhauer M., et al. 2007, MNRAS, submitted 

Fukugita M., Ichikawa T., Gunn J. E., Doi M., Shimasaku K., 

Schneider D. P. 1996, AJ, 111, 1748 
Koposov, S. E., Glush kova, E. V & Zolotukhin, I. Y. 2007, ArXiv 

e-prints, 709, arXiv:0709 . 12751 accepted to A&A 
Gunn, J.E. et al. 1998, AJ, 116, 3040 
Gunn, J.E. et al. 2006, AJ, 131, 2332 

Gnedin, O. Y., Zhao, H., Pringle, J. E., Fall, S. M., Livio, M., & 

Meylan, G. 2002, ApJ, 568, L23 
Harris W.E., 1996, AJ, 112 1487 

Huxor, A. P., Tanvir, N. R., Irwin, M. J., Ibata, R., Collett, J. L., 
Ferguson, A. M. N., Bridges, T., & Lewis, G. F. 2005, MNRAS, 
360, 1007 

Hogg, D.W., Finkbeiner, D.P., Schlegel, D.J., Gunn, J.E. 2001, AJ, 
122, 2129 

Irwin M. J., Hatzidimitriou D. 1995, MNRAS, 277, 1354 
Irwin M. J., et al. 2007, ApJ, 656, L13 

Kleyna, J., Wilkinson, M. I., Evans, N. W., Gilmore, G., & Frayn, 

C. 2002, MNRAS, 330, 792 
Klypin A., Kravtsov A. V., Valenzuela O., Prada F. 1999, ApJ, 

522, 82 



Klypin, A., Zhao, H., & Somerville, R. S. 2002, ApJ, 573, 597 
Koposov, S., & Bartunov, O. 2006, Astronomical Data Analysis 

Software and Systems XV, 351, 735 
Koposov, S., Bartunov, O., Belinskiy, A., Karpov, S. 2007, 

Astronomical Data Analysis Software and Systems XVI, 376, 

34 

Koposov, S., et al. 2007, ApJ, 669, 337 

Lindenberg, T., 1993, International Journal of Computer Vision, 
11 , 3, 283 

Lindenberg, T., 1998, International Journal of Computer Vision, 
30, 2, 79 

Lupton, R., Gunn, J., & Szalay, A. 1999, AJ, 118, 1406 

Lupton, R., Gunn, J. E., Ivezic, Z., Knapp, G. R., & Kent, S. 2001, 

Astronomical Data Analysis Software and Systems X, 238, 269 
Mackey, A. D., & Gilmore, G. F. 2004, MNRAS, 355, 504 
Martin, N. F., Ibata, R. A., Chapman, S. C, Irwin, M., & Lewis, 

G. F. 2007, MNRAS, 380, 281 
Metz, M., Kroupa, P., & Jerjen, H. 2007, MNRAS, 374, 1125 
Moore B., Ghigna S., Governato F., Lake G., Quinn T., Stadel J., 

Tozzi P. 1999, ApJ, 524, L19 
Moore B., Diemand J., Madau P., Zemp M., Stadel J. 2006, 

MNRAS, 368, 563 
Odenkirchen M., et al. 2001, ApJ, 548, L165 
Page, L., et al. 2007, ApJS, 170, 335 

Pier, J.R., Munn, J. A., Hindsley, R.B, Hennessy, G.S., Kent, S.M., 

Lupton, R.H., Ivezic, Z. 2003, AJ, 125, 1559 
Piotto, G., et al. 2007, ApJ, 661, L53 

Schlegel, D. J., Finkbeiner, D. P., k, Davis, M. 1998, ApJ, 500, 
525 

Silverman, B. W. 1986, Monographs on Statistics and Applied 

Probability, London: Chapman and Hall, 1986, 
Smith, J. A., et al. 2002, AJ, 123, 2121 
Somerville, R. S. 2002, ApJ, 572, L23 

Stoehr F., White S. D. M., Tormen G., Springel V. 2002, MNRAS, 
335, L84 

Stoughton, C. et al. 2002, AJ, 123, 485 



14 



Koposov et al. 



Tucker D., et al. 2006, AN, in press 

de Vaucoulcurs, G., de Vaucouleurs, A., Corwin, H. G., Jr., Buta, 
R. J., Paturel, G., & Fouque, P. 1991, Volume 1-3, XII, 2069 
pp. 7 figs.. Springer- Verlag Berlin Heidelberg New York, 

van den Bergh, S. 2000, The galaxies of the Local Group, by 
Sidney Van den Bergh. Published by Cambridge, UK: Cambridge 
University Press, 2000 Cambridge Astrophysics Series Series, vol 
no: 35, ISBN: 0521651816., 

Walsh, S. M., Jerjen, H., & Willman, B. 2007, ApJ, 662, L83 

Whiting, A. B., Hau, G. K. T., Irwin, M., & Verdugo, M. 2007, 
AJ, 133, 715 

Willman, B., Dalcanton J., Ivezic Z., Jackson T., Lupton R., 
Brinkmann J., Hennessy G., Hindsley R. 2002, AJ, 123, 848 



Willman, B., et al. 2005a, ApJ, 626, L85 
Willman, B., et al. 2005b, AJ, 129, 2692 
York D. G., et al. 2000, AJ, 120, 1579 

Zinncckcr, H., Keablc, C. J., Dunlop, J. S., Cannon, R. D., 
& Griffiths, W. K. 1988, The Harlow-Shapley Symposium on 
Globular Cluster Systems in Galaxies, 126, 603 

Zucker D. B., et al. 2006a, ApJ, 643, L103 

Zucker D. B., et al. 2006b, ApJ, 650, L41 



